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Abstract 

We consider a wide class of model equations, able to describe wave propagation in 
dispersive nonlinear media. The Korteweg-de Vries (KdV) equation is derived in this 
general frame under some conditions, the physical meanings of which are clarified. 
It is obtained as usual at leading order in some multiscale expansion. The higher 
order terms in this expansion are studied making use of a multi-time formalism and 
imposing the condition that the main term satisfies the whole KdV hierarchy. The 
evolution of the higher order terms with respect to the higher order time variables can 
be described through the introduction of a linearized KdV hierarchy. This allows one 
to give an expression of the higher order time derivatives that appear in the right hand 
member of the perturbative expansion equations, to show that overall the higher order 
terms do not produce any secularity and to prove that the formal expansion contains 
only bounded terms. 

1 Introduction 

Soliton theory is based on the resolution of the so-called "completely integrable equations" 
by means of the inverse scattering transform (1ST) method. When they describe particular 
physical situations, these equations arise as asymptotics of well-established models, and are 
derived from the latter by means of a multiscale expansion, or some equivalent formalism. 
The integrable nonlinear evolution equations appear as the leading order approximation 
in this perturbative approach. Corrections to this first approximation are often to be 
considered. At a given propagation time (or distance), more accuracy can obviously be 
obtained by retaining more than the main term only in the power series. Such additional 
terms are the "higher order terms". On the other hand, for a very long propagation 
time, the physical solution goes away from the theoretical soliton because the nonlinear 
evolution equation gives only a first order approximation of this evolution. Corrections to 
the equation must be taken into account. The partial differential equations giving such 
corrections are what we call the "higher order" ones. 

Consideration of such corrections have given rise to perturbation theories for solitons |Hj 
that have found applications, e.g. in the frame of the physics of optical solitons in fibres, 
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which are related to optical telecommunications @]. The effect of the higher order equa- 
tions appears in these theories notably through the renormalization of the soliton velocity. 
From the purely theoretical point of view, it has been recently shown that the solvability 
of the higher order equations can in some sense be related to the problem of the complete 
integrability of the basic system 110] . Furthermore, a physical interpretation of the 
equations of the Korteweg-de Vries (KdV) hierarchy has been found by Kraenkel, Manna 
and Pereira in the framework of the theory of water waves 1121 113] . 

The Korteweg-de Vries (KdV) equation, which is known to describe the propagation of 
long waves in shallow water [2] , is the first equation that was solved by the 1ST method |S[ 
14;. More physically speaking, the first observed soliton was a long wave in shallow water. 
The KdV equation appears as the evolution equation satisfied by the dominant term of the 
quantity describing the wave in some multiscale expansion. We investigate in this paper 
the higher orders terms of such an expansion. Our results apply a priori to the Maxwell- 
Landau model that describes electromagnetic wave propagation in ferromagnetic media, 
but, because they do not involve the explicit computations particular to this situation 
(written down elsewhere |16j). they are presented in a general abstract frame. The results 
are thus expected to apply to other physical situations. The conditions under which the 
KdV model can be derived in this general frame are discussed below from both the physical 
and the mathematical points of view. Regarding the higher order equations such conditions 
cannot be written down explicitly in a completely satisfactory way from the mathematical 
point of view, but their physical meaning is clarified. Further it has been proved ^0] 
that these conditions are satisfied by the Maxwell-Landau model which describes wave 
propagation in ferromagnets. 

We use, as did Kraenkel et al, a multiple time formalism. These authors have shown 
that the evolution of the dominant term relative to the higher order time variables is 
governed by the KdV hierarchy. Kodama and Taniuti showed that the evolution of the 
higher order terms relative to the first time variable is governed by the linearized KdV 
equation, but the way in which the higher order terms depend on the higher order time 
variables was not yet clarified and it was not proven that the higher order terms do not 
produce secular terms, while the coherence of the multiscale expansion necessitates the 
elimination of any secularities. We find that the evolution of these terms can be described 
by means of some linearized KdV hierarchy. This is the main result of the paper and 
enables us to prove that no unbounded term does appear in the expansion. 

The paper is organized as follows. In Section 2 we describe the frame of the multiscale 
expansion, derive the KdV equation and discuss the physical hypotheses that allow this 
derivation. In Section 3 the problem of secularities is presented and the introduction of the 
KdV hierarchy allows us both to remove the secularity producing terms due to the main 
order and to determine the evolution of the latter relative to the higher order variables. 
Section 4 is devoted to the higher order time evolution of the higher order terms and to 
the linearized KdV hierarchy. Section 5 contains a conclusion. The linear part of the right 
hand member of the equations of the perturbative expansion plays an important role in 
the treatment of the secularity producing terms. Therefore the corresponding coefficients 
are computed explicitly in an appendix. 



Higher Order Terms in Multiscale Expansions: A Linearized KdV Hierarchy 



327 



2 A multiscale expansion 



2.1 The model and the scaling 



We consider some set of partial differential equations that can be written as 



(d t + Ad x + E) u = B(u,u) 



(2.1) 



where the function u of the variables x and t is valued in M. p , A and E are some p x p 
matrices, and B : W x W — > W is bilinear. This system can describe the propagation of 
electromagnetic waves in a ferromagnetic medium, according to the Maxwell and Landau 
equations [El- We consider the abstract system (|2.1jl rather than this latter particular 
situation for sake of simplicity. Indeed the study of a specific case would imply explicit 
computation of all the quantities involved by the multiscale expansion and increase consi- 
derably the size of the expressions, while the matter of the present paper does not depend 
upon the particular form of these quantities. It is in fact rather general, according to 
the fact that KdV solitons arise in many other physical frames. Considering an abstract 
frame thus avoids much computational detail, but in fact is not so easy because the deriva- 
tion of the KdV model is not always achievable. It involves several assumptions, some 
of which are rather strong. These assumptions are introduced at the points where they 
have to be used. Then in subsection 12.31 we give their physical interpretation. In any case 
the required assumptions are satisfied by the Maxwell-Landau model, describing waves 
in ferromagnets. The proof of this is given in [IB] , Thus the results of the paper are 
valid and completely justified for this latter model. The use of the abstract system (|2.1|) 
presents the two following advantages. Firstly it simplifies the algebra because it avoids 
considering many details peculiar to the physics of ferromagnets. Secondly it allows a dis- 
cussion of the physical conditions required for the occurrence of KdV solitons, despite 
some feature concerning the higher order equations that cannot be completely solved from 
the mathematical point of view. 

Regarding the system (|2,lj) the following assumptions are usually made: A is assumed 
either to be completely hyperbolic or symmetric and E is assumed to be skewsymmetric. 
The two latter hypotheses are closely related to the conservative character of the system 
and are satisfied by the Maxwell-Landau system under the scalar product 



(with the notations of ^H]). Note that this scalar product depends upon the zero-order 
term. Only weaker assumptions are necessary for the formal derivation of the KdV equa- 
tion, as we will see below. However, some essential symmetry property of the higher order 
equations is satisfied only when the system is conservative. This is ensured, regarding the 
linear part of the system, by the symmetry and skewsymmetry hypotheses above. The 
corresponding assumption on the nonlinear part cannot be expressed in a simple way. 
The variable u is expanded in a power series of some small parameter e as 



( 




E 2 + H 2 + olM 2 



u = e 2 u 2 + e 3 u 3 H 



(2.2) 



and the slow variables 



£ = £ {x - Vt), t = e 3 t 



(2.3) 
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are introduced. V is a speed to be determined. This multiscale expansion is well-known. 
It yields the Korteweg-de Vries (KdV) equation at leading order, e 2 , in many physical 
cases. Starting the expansion at order e 2 corresponds physically to specifying a low value 
for the order of magnitude of the amplitude. Setting a term of order e in expansion (|2.2|) 
is a priori possible, but it is easily checked that it does not yield any evolution equation 
regarding the chosen time scales. Therefore we can omit it. It is seen below that, as 
is described in [HI El, secularities, i.e. linear growth of the solutions with time, appear 
in the higher order terms. These secularities must be removed and this is achieved by 
imposing to the leading term some particular dependency with regard to higher order 
time variables (the equations of the KdV hierarchy) . Therefore we introduce the variables 
T2, T3, . . . which are defined by 

Tj = e^ +1 t (j > 1). (2.4) 

With this notation r = t\. 

2.2 The Korteweg-de Vries equation 

At order e 2 equation (|2.1|) yields 

Eu 2 = 0. (2.5) 

Thus U2 belongs to the kernel, ker(E), of E. We assume that the range, Hg(E), of E is 
in direct sum with its kernel (hypothesis 1). This is satisfied if E is skewsymmetric. We 
also denote by Ilo and Qo respectively the projectors onto kei(E) parallel to Rg(E) and 
onto Jig(E) parallel to ker(£'). Equation (|2,5|) can be written as 

u 2 = n U2- (2.6) 

At the following order, e 3 , we get the equation 

Eu 3 + (A - V)d/:U2 = 0. (2.7) 

The projection of equation (|2.7|) on ker(E) is 

n (A- V)d^U2 = 0. (2.8) 

Because of (|2.fi[) . equation (|2.8|) has a nonzero solution if V is an eigenvalue of Ilo Alio. 
We assume that V is a simple nonzero eigenvalue of Ilo Alio, an d call ao an eigenvec- 
tor, 111 and Qi the associated projectors: IT the projector on aoM, Q\ the projector on 
Rg (Ilo (A — V")Ho), so that IT + Q\ = Hq. This definition assumes that W is the direct 
sum of the range and the kernel of IIo(A — V)IIo. This implies that not only the eigenspace 
but the characteristic space of the operator iloArio relative to the eigenvalue V has di- 
mension 1 (hypothesis 2). The assumption of the complete hyperbolicity of A ensures 
that this hypothesis is satisfied. When A is assumed to be symmetric, the dimensions of 
characteristic and eigenspace are the same. 

According to equation (|2.8j) d^U2 G ker(IIo(A — V)IIo). Thus 112 = ao^2, where the 
function if 2 has to be determined. The Qo-P r ojection of equation 1)2.7(1 is 



Q u 3 = aid^(p 2 . 



(2.9) 
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The vector coefficient a\ is defined by a\ = — E~ l (A — V)clq, where E~ l is a partial inverse 
of E. E~ 1 is precisely defined as follows. We call E the restriction and corestriction of 
QoEQo to Kg(E): E is invertible. We denote by the zero operator defined on ker(E). 
Then E^ 1 is the direct sum of and E -1 . In other words, the part of the matrices of E 
and E~* corresponding to Kg(E) are inverses one of each other and the matrix of E~ l is 
completed by zeros to make a p x p matrix. 
The equation of order e 4 is 

Eu 4 + (A - V)d^u 3 = B(u 2 ,u 2 ). (2.10) 

It is divided into 3 parts by using the projectors Qo, Q±, IIi. The Hi-projection is a com- 
patibility condition for (p 2 . Using IIi (^4 — V)Ilo = 0, and splitting U3 into U3 = H0U3+Q0U3, 
we find that the projection reduces to 

qd}?2 = rtpl (2.11) 

with q = Tl±(A — V)a\ and r = IliB(ao, ao). If (q,r) ^ (0,0), equation (|2.11j) has no 
bounded nonzero solution. The following conditions must thus be satisfied (hypotheses 3 
and 4), viz 

niS(a ,a ) = 0, and Hx(A - V)a x = 0. (2.12) 

The Q\- and then the Qo-P r ojection of equation (|2.10j) yield expressions for Q1U3 and Q0U4, 
respectively, as 

Qiu 3 = a' x d^2 + ai / f2 2 , (2.13) 

Q0U4 = aid^tpz + a 2 d^(p 2 + a 2 f2 2 - (2-14) 

Denoting by (A — V)^ 1 a partial inverse of Ho(j4 — U)no defined in an way analogous 
to E , we define the vector coefficients by 

a' 1 = -(A-V)- 1 Q 1 (A-V)a 1 , (2.15) 

a 1 = (A-V)- 1 Q 1 B(a ,a ), (2.16) 

a 2 = -E- 1 Q (A-V)(a 1 + a[) , (2.17) 

OL2 = E- X Q Q (B(oo, ao) - (A - V>i) . (2.18) 

With the use of all previous results the compatibility condition (ni-projection) of the 
equation of order e 5 yields the following evolution equation for ip 2 : 

d T (f 2 + 0<P2d£<p 2 + 7<9f (f2 = 5tf2 i f2 2 - (2.19) 

The scalar coefficients (5, 7, 5 are defined by 

a (3 = 2IIi(A - V)oL 2 - 2U 1 B(a , a x + a[), (2.20) 
a j = n 1 (A-V)a 2 , (2.21) 
a 5 = 2Il 1 B(ao,ai). (2.22) 
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If the additional condition (hypothesis 5) 

11x5(00, ai) = (2.23) 
is satisfied, equation Q2.19[) reduces to the KdV equation 

d T ip 2 + f3ip 2 d^ 2 + 7df ¥>2 = 0. (2.24) 

2.3 Physical meaning of the hypotheses 

The above formal derivation of the KdV equation necessitates 5 hypotheses that have been 
written down explicitly. Two of them, hypotheses 3 and 5, involve the nonlinearity, the 
three other the derivatives only. A physical interpretation of the latter can be found in 
the dispersion relation of the linearized system 

(d t + Ad x + E)u = 0. (2.25) 

The pulsation and polarization vector corresponding to a given wave vector, k, are denoted 
by uj(k) and u{k) respectively. They satisfy 

(-ioj(k) + Aik + E) u(k) = 0. (2.26) 

The long wave approximation corresponds to k = and w(0) = 0. Then equation (|2.2(i|) 
reduces to Eu(0) = 0. It is equation (|2,5|) . Taking the derivative of equation (|2.26l) with 
regard to k and then setting k = we obtain 

i[A - o/(0)]u(0) + E u'(0) = 0. (2.27) 

The prime denotes the derivative with regard to k. Together with the interpretation of 
the long wave approximation as a limit of oscillating waves for k close to 0, the analogy 
between equation (|2.27|) and equation (|2.7|) allows the identification between u(0) and ao 
on one hand, and between u/(0) and V on the other. V appears as a long- wave limit of 
the group velocity. Hypothesis 2 is thus that a unique polarization can propagate with 
this velocity. If this assumption is not satisfied, interactions between the various waves 
with same velocity must be taken into account. 
Equation IJ2.27J) yields also 

Q o u'(0) = -E- l i[A - V]a . (2.28) 

Thus Qqu'(0) = ia%, with the definition above of a%. Taking once again the derivative of 



equation (|2.26j) with regard to k we obtain, for k = 0, 

-iuj"(0)a + 2i[A - V]u'(0) +Eu"(0) = 0. (2.29) 

Taking the Ili-projection of equation (|2.29j) . decomposing u'(0) into u'(0) = ia% + Hou'(0), 
and taking into account the relation ni(A — y)Ilo = 0, we get 

u/'(0)a = 2ffli [A - V]m. (2.30) 
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Hypothesis 5, in (|2.12jl . is thus 

d 2 uo , , 

The following feature can be observed: u/"(0) can be computed by the same method, using 
the second derivative of equation 1)2. 26|) . The expression obtained gives an interpretation 
of the coefficient 7 of the third derivative in the KdV equation, 1)2.19)1 . as 

— 1 dPuj , . , 

Hypothesis 1 appears mainly as a technical assumption without special physical mean- 
ing needed to solve equation 1)2.28)1 . If E is skewsymmetric, this hypothesis is satisfied. 
Assume that E has a nonzero symmetrical part E s . Taking the scalar product of equation 
()2.26)1 by u(k), then comparing the result to its complex conjugate and making use of the 
symmetry assumptions (A being symmetric), we obtain the following expression for the 
imaginary part uJiik) of u(k): oJi{k) = —{u\E s u)/{u\u). Thus the system cannot be con- 
servative if E is not skew-symmetric, and hypothesis 1 is always satisfied by conservative 
systems. 

The two hypotheses 3 and 5 concerning the nonlinearity obviously cannot be understood 
through the study of the linearized equation 1)2.25)1 . Hypothesis 3 can be written as 

niB(ni,ni) = o, (2.32) 

while the assumption 

ni5(IIi,Qi) =0 (2.33) 

implies hypothesis 5. Conditions 1)2.32)1 and 1)2.33)1 are particular expressions of a very 
general condition called the "transparency" condition in the rigorous mathematical the- 
ory of multiscale expansions by Joly, Metivier and Rausch . Condition (j2.32j) excludes 
quadratic self-interaction for the chosen propagation mode, while condition ()2.33)1 ex- 



cludes interaction at the same order for different polarizations. If condition 1)2.32)1 . that 
is hypothesis 3, is not satisfied, the nonlinear term appears sooner in the expansion. The 
nonlinear evolution of the wave will be described by a nonlinear evolution equation other 
than the KdV equation and for other space, amplitude and time scales. If condition 1)2.33)1 . 
or more precisely hypothesis 5, is not satisfied, the above computation is valid, but the 
evolution equation is equation 1)2.19)1 instead of KdV. It is an integro-differential equation 
involving a cubic nonlinear term. 



2.4 The linearized KdV equation 

Regarding the higher order terms, the evolution equation relative to the first order time 
r = n is found in the following way. Equation 1)2.1)1 . written for some given order e n , 
yields a compatibility condition, its Hi-projection 

n-2 

III (A - V)d i Q u n ^ 1 +^2a rj U 1 u n - 2 j-i = Y, U i B ( u P' u n- P ), ( 2 -34) 
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and recurrence formulas for the determined parts of u n , its Qo-P r °j ec ti° n 

n-2 



Qou n — E 1 Qq 
and its Qi-projection 



-{A-V)d i u n . l -Y,dr ] 



u 



n-2j-l + 

p=2 



j n—p ) 



(2.35) 



Q lUn = (A-V)- L Q 



(A — V)QoU n — / 5 Tj n„_2i + ^ / B{ 

j^l ^ p=2 ^ 



(2.36) 



The component of u n belonging to Rg(Ili) is proportional to some unknown real func- 
tion cp n , and u n is 



<Wn + QoUn + QlU r , 



(2.37) 



We call S p the set of expressions involving 112,113, ... ,Up, but not u p+ \ and subsequent 
terms. The same notation S p will hold below for any function belonging to this set. Due 
to equations Q2.35J) and Q2.36JI u n can be written as 



<Wn + S. 



n— 1 • 



(2.38) 



The functions ip n are the unknowns of the problem. Note that, with these notations, S p is 
the set of expressions involving cp%, ^3, . . . , cp p and their derivatives. Applying two times 
the recurrence formulas (|2.35|) and (|2.36|) onto the expression (|2.38|) of u n we obtain 



O-Ofn + QoUn + QlU n , 



with 



Q U n = aid^(p n -i + a 2 9|(/?n-2 + ^a 2 ip2<fn-2 + Sn- 

Qiu n = a! x d^ip n ^i + 2a\ j <p2Vn-i + <S n -2- 



(2.39) 

(2.40) 
(2.41) 



Equations (|2,39j) - (|2,41j) are correct if n — 2 > 2, i.e. n ^ 5. For n = 4, the coefficient 2 
before a 2 vanishes, as seen above. The dependency of Q\u n with regard to tp n -2 is not 
needed in the computation of the higher order equations below. This is fortunate because 
of several integral nonlinear terms arising in it that would greatly hinder the computa- 
tion. When use is made of the expressions (|2.39|) to (|2.41|) of u n in the compatibility 
condition (|2.34j) . it gives an evolution equation for <p n -3. Indeed, u n does not appear in 
the equation. The dependency with regard to (p n -i, tfn-2 and ip n -z is explicitly written. 
Because of the condition IIi(.A — V)a = 0, which determines the velocity V, <p n —i vanishes 
from the equation. Because of the relations (|2.12|) and (|2.23l) that are satisfied under the 

present hypotheses, (f n -2 and a term 5if2 ( / ¥>2¥>n-3 + / V2 2 J also vanish and the 



(2.42) 



equation obtained is, with n — 3 = I: 

d T1 ipi + (3dj: (ip2Vi) + jdlipi = Ei (ip 2 , ipi-i) . 
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Equation (|2.42j) is the linearized KdV equation obtained by the linearization of the KdV 
equation (|2.24l) about its solution if 2 with an additional right hand member H;(</?2,..., 
f l— 1) depending on the previously determined terms. 

This is true with no further hypothesis than the existence of some simple eigenvalue V 
of no (A — V)Uq and conditions (j2.12j) and (|2.23|) . However, in the general case, due to 
the term B(u p , u n +i- p ) in the recurrence formula (|2.36|) . the right hand member H/ 
involves a priori many integrations relative to the variable £. The property 



satisfied in ferromagnets, ensures their vanishing. 

3 The second order time evolution and the KdV hierarchy 
3.1 The problem of the secularities 

It is well known that the KdV equation ()2.24j) is completely integrable, i.e. that the 
Cauchy problem for it can be solved by use of the Inverse Scattering Transform (1ST) 
method The 1ST method gives also explicit formulas for the resolution of the 

linearized KdV equation (|2.42|) P^llBj. This latter equation is linear, but it is necessary to 
solve it by this method, because the solution if 2 of the KdV equation intervenes in it as an 
essential parameter. In the general case <f2 can be only expressed in terms of its inverse 
transform, thus also the solution of the linearized KdV equation (J2.42|) . Unfortunately 
the solution computed this way is not always bounded (I mean uniformly bounded in £ 
as r — > +00). As an example, assume that at r = 0, tpi = 0, and replace the right hand 
member (<p%, . . . , <£z-l) by d^f2- Then the solution ifi of equation (|2.42|) is [HI E] 



A solution like (|3.1j) is called "secular" . This phenomenon occurs when a term in the right 
hand member "resonates", that is, from the mathematical point of view, is a solution of 
the homogeneous equation ( (|2.42|) with S; (932, • • • , fi-i) = 0). Physically, this resonance 
phenomenon occurs when each component of the inverse scattering transform of the source 
term (some part of the right hand member {f2, <£3, ■ ■ ■ > tpi—l) ) evolves in time in the 
same way as the corresponding component in the transform of the main term. In order 
to remove the secular terms, we use the method of Kraenkel, Manna, and Pereira |121 113j 
that consists in the introduction of additional evolution equations, which describe the 

evolution of the main term relative to the higher order time variables T2, T3, We 

check that, according to the papers cited, these equations must be those of the so-called 
KdV hierarchy, normalized in order to cancel the linear terms in the right hand member 
3; {if 2, f 3, ■ ■ ■ , fi-i)- Indeed references (HI El state that the secularity producing terms 
are the linear ones. We showed in |15j that they are rather the derivatives of the conserved 
densities of the KdV equation, and that the procedure of Kraenkel et al. in fact removes 
completely these latter terms from the right hand member. 

3.2 The next order equation 

The standard expansion used to derive the KdV equation in hydrodynamics has some 
parity relative to e. As an example, Su |17j uses series expansions in integer powers of 



QlB{;-)=0, 



(2.43) 



fi = rd^f2- 



(3.1) 
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some perturbative parameter e, defined in such a way that e = e 2 , e being the perturbative 
parameter of the present paper. This corresponds to a vanishing of all odd terms in the 
present expansion. The equation (|2.42|) governing the t\ evolution of (p n is obtained at 
order e n+3 , which is even when n is odd. Due to the homogeneity properties of the 
expansion, the right hand member 3 n (992, ■ ■ ■ , l -p n -i) contains only even order derivatives 
of </?2i which describe losses. Physically, if the initial system ()2.1|) is conservative, the 
same feature can be expected for any equation in the expansion. Therefore the even order 
derivatives of 992 must vanish. We assume that the present model has the same property 
in the following sense: H3 (932) = 0. Thus, with a zero initial condition, 993 = 0, and so 
on: if all (p p with odd p are zero up to (n — 2), n being odd, then E n (992, ■ ■ ■ , tpn-i) is 
zero and, with a zero initial condition, ip n is also zero. Thus, if the corresponding initial 
conditions are zero, all the H n and all the (p n with odd values of the integer n are zero. 
This occurs in the particular case in which the system (|2,1|) describes wave propagation 
in a ferromagnetic medium. This last feature is proved, using homogeneity properties, 
in [T^]. Few results can be proved on this point in the general case. The right hand 
member H3 (992) contains a term Kdf(p%, where the coefficient K is given by 

a K = U^A - V)E~ 1 Q (A - V) (l - (A - V)~ X Q\(A - V)) a 2 . (3.2) 

In the same way as the expression (|2.3()|) for uj" has been derived, we can give an expres- 
sion for a/ 4 ) = d 4 uj / dk 4: \ k _Q. After comparison with 1)3.2)) . it is seen that u/ 4 ) = 24iK. 
Observe that the right hand member of 1)3. 2 j) contains only real terms. Thus K is real and 
u/ 4 ) is purely imaginary. The requirement that the linear term Kdf(f2 in H3 (992) vanish 
is equivalent to the requirement that u/ 4 ) be real, which must be satisfied if the system is 
conservative. This generalizes to the higher even order derivatives and ensures the vanish- 
ing of the linear terms, but we omit the proof because we must admit the generalization 
of this feature to the nonlinear terms. Note that only the linear terms are suspected to 
be secularity producing. Thus the vanishing of the linear term is expected to ensure that 
all terms in the expansion are bounded. Nonetheless, we make the assumption above, 
because it is satisfied in all the examples of which we know, even if we are not able to 
prove it in the general case. 

With this hypothesis the second nontrivial equation of the perturbative expansion is 
the equation p. 42)) for 994. Its right hand member H4 (992) is polynomial with regard to 992 
and its derivatives, with the homogeneity of the terms of order e 7 in the expansion and 
can thus be expanded according to formula ()A.16)) as 

^(^2)= s (K-);>i>*0lI / ^ %>2 + o 2 . (3.3) 

i>1 / 

Ej^i 2jmj+k=5 

{O2 designates here an expression in 992 without linear terms.) The indices (jrij)j^i, k in 
this sum can take the values 

(K)j>i,A;) = ((0,...),5); ((1, 0, . . .), 3); ((2, 0, . . .), 1); ((0, 1, 0, . . .), 1) (3.4) 

only. Furthermore, the coefficient H ((0, 1,0, . . .), 1) of d T2 (f2 is —1. Thus expression 1)3.3)) 
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can be expanded in the following way: 



s 4 (¥>2) = -a^2 + 3((2,o,...),i) .O dm 

+ H ((1,0, . . .), 3) (J d T1 \ df V2 + 3 ((0, . . .), 5) dl<p 2 + 2 . 



(3.5) 



As (f2 satisfies the KdV equation (|2.24j) . we have 

/ ^^ = -7^ + 02. (3.6) 

J — oo 

Thus 

^4 (<P2) = -dn<P2 ~ 72«9|</? 2 + 2 , (3.7) 

with 

-72 = S ((2, 0, . . .), 1) (- 7 ) 2 + S ((1, 0, . . .), 3) (-7) + S ((0, . . .), 5) . (3.8) 
The coefficient 72 is explicitly computed in the Appendix (equation (|A.19|) ). 

3.3 The KdV hierarchy 

The linear terms computed above are secularity producing. Because the solution 934 of 
the linearized KdV equation must be bounded, they must vanish. This imposes some 
evolution equation for ip 2 relative to the second-order time variable r 2 , such that 

3 T2 9?2 = -723fv?2 + 02. (3.9) 

The nonlinear terms in equation (|3.9j) are not free, but imposed by the compatibility con- 
ditions between the KdV equation (|2.24|) and (|3.9j) : the Schwartz conditions dr\dT 2 ^ 2 = 
dT 2 dr\if> 2 . Kraenkel, Manna, and Pereira |13j have shown that the only equation that has 
the same homogeneity properties as H4 and that satisfies this condition is the second equa- 
tion of the so-called KdV hierarchy. The same requirements are found at higher orders. 
The KdV hierarchy is the following set of equations [2] : 

&T n v = dxC n v (n integer), (3.10) 

where 

£ = - v + \f dX(d x v). (3.11) 

For n = 1 it is the KdV equation, but with values of the coefficients (3 = 3/2, 7 = 1/4. 
The coefficients become equal to those of equation (|2.24(l by setting 

v = ^-ip 2 , X = i and Ti=4 7 r 1 . (3.12) 

67 
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For n = 2 the equation Q3.10JI of the hierarchy is 

9t 2 v = ^d x v + \{d x v)d 2 x v + h -vd x v + ^v 2 d x v. (3.13) 

The important feature is the existence of the Hirota r- function [2], that is a function of 
all the variables (X, Ti, T 2 , . . .), related to v by 

v(X,T 1 ,T 2 ,...) = 2d x \nr(X,T 1 ,T 2 ,...) (3.14) 

(take care to avoid confusion between the Hirota r-function and the time variables Tj ) . The 
existence of r insures that a solution v of the whole system yielded by all equations of the 
hierarchy exists. As the system admits a solution, the Schwartz conditions are satisfied 
at any order, i.e. d^d^y = d^d^v for any j, p ^ 1. Furthermore these conditions 
are satisfied formally and identically by the equations themselves and not only for some 
particular solution. The authors cited have checked by formal computation on the first 
orders that, apart from the symmetries of the KdV hierarchy, it is the only compatible 
system which has the required homogeneity properties and for which the first equation 
is the KdV one. The symmetries of the hierarchy are, after those of the KdV equation, 
a free scaling coefficient for each time variable. We must identify 

T 2 = -I672T2 (3.15) 

with 72 given by equation (|3.8|) and impose that ip2 satisfy the evolution equation 
-1 „ 

d T2 if2 = d^C ip 2 (3.16) 



-1*2 P.. , P f 6 



I672 

with 

c --&--fr* + £;L' mfP,) - (3J7) 

Then H4 {^2) n ° longer contains any linear term. This implies, due to the procedure, that 
it no longer contains any secularity producing term |15j . 

The linear terms in </? 2 are removed in the same way at each order. We impose for each 
p ^ 2 that 

P^ T ^ 2 = d^ 2 , (3.18) 
(C as above) with 7 p defined by 71 = 7 and 

p-i 

7 P +i= E s (K)i^-i> fc )II(-^) mi - ( 3 - 19 ) 

7=1 

£P- 1 1 2jm,+fc=2p+3 

For p = 1 equation (|3.18|) coincides with equation ()2.24j) with j3\ = (5. We use here the 
same scheme as in the first case p = 2. We have 

CP=(^d 2 X + O l . (3.20) 
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Equation ()3.18j) gives 

d Tp P2 = -l P df +1 y 2 + 02. (3.21) 

Then the linear term in <p2 vanishes from 3 p (</?2, <pi, . . . , </?p-2)- I n this way all secularity 
producing terms due to ip2 vanish at any order. 

4 The third order and the linearized KdV hierarchy 
4.1 The higher order time evolution of the higher order terms 



The third equation of our perturbation expansion is equation 1)2. 42|) . written for 1 = 6, 
with its right hand member E.q ((/?2,<£>4) given by formula ()A,16)) . It is of order e 9 . It 
involves not only tp2, tp4 and their space derivatives but also d T1 (p2, d T2 if2, d T3 ip2, d Tl ip4 
and d T2 tpi of orders e 5 , e 7 , e 9 , e 1 and e 9 respectively. d T1 if2 is replaced by space derivatives 
and nonlinear terms using the KdV equation 1)2 .24|) . We make use of the second and third 
equations of the hierarchy for cp2 so that d T2 (f2 and d Ti <p2 also vanish, replaced by functions 
of the space derivatives. The right hand member becomes 

^5(^2,^4) = <9 T2 c/J 4 + ^((2,0, ...),1) (^j d T1 ^j 8^4 

+ S ((1, 0, . . .), 3) d Tl \ dl<p 4 + 3 ((0, . . .), 5) 9^4 + 2 . (4.1) 

(Here 02 represents some expression in <p2, <P4, their derivatives and primitives, without 
linear terms). The ri-evolution of ip4 is also known, and described by the linearized KdV 
equation 1)2. 42|) . with the right hand member 34(^2)- Using this relation removes the 
explicit dependency on d T1 (p4, and the right hand member (if2, 924) reduces to 

3e (v>2, ^4) = 9 T2 (f4 + 72^1^4 + 02- (4.2) 
Two questions arise at this point: 

• How is the T2-dependency of 924 defined? 

• Is the right hand member (^2, ^4) secularity producing due to the terms in 934, 
in particular the linear ones? 

if 4 is the solution of the linearized KdV equation 1)2.42)1 . with the right hand mem- 
ber 34 ((^2) and some given initial data (p 4 (£, t\ = 0). It is expressed as an integral and 
linear combination of the squared Jost functions, defined in the solution of the KdV equa- 
tion through the 1ST method [HJ E] • We denote by cp^ the solution of the homogeneous 

linearized KdV equation with the same initial data V?4 (£, t\ = 0) = <p^, t\ = 0), and 

(2) — 
by ^4 the solution of the linearized KdV equation with the right hand member 34 yf2) 

and vanishing initial data <p^ (£, t% = 0) = 0, so that 994 = ip^ + (p? . It is shown in 15_ 

that (y?^ is secularity producing while ip± is not. 

Since the T2-dependency of the solution ip2 of KdV is determined, the squared Jost 

functions are also known for all (£, n, T2), and so are the spectral components of the right 
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^ (2) 

hand member H4 [if 2) ■ Thus <p\ is completely determined and its T2-dependency is well 
defined without additional condition. Furthermore the evolution equation that makes 
explicit this dependency is not needed for the computation of the function y?4(£,Ti,T2). 
The situation is different with ip^ , while the initial condition is a priori free. In the 
previous subsection we saw that the compatibility conditions fixed the T2-dependency 
of if2, apart from a scaling coefficient for this time variable. The T2-dependency of ipr± is 
determined below in the same way, but the scaling coefficients are no longer free. 

Consider some solution vq of the KdV equation, and v\ of the homogeneous linearized 
KdV equation. We use here the normalization of formulas (|3.10j) - ()3.11j) . The latter 
equation is 

dr^i + ^dxivoVi) + ^d xVl = 0. (4.3) 

We assume also that vq and v\ are smooth functions of all the variables X, T\ , T2, . . ., and 
that their dependency relative to each time variable satisfies the homogeneity properties 
of the KdV hierarchy. We saw above that under these conditions the function vq, solution 
of the KdV equation, satisfies the complete hierarchy, with some scaling constant for each 
time variable. Consider some small parameter ij. The function 

v = vo + ijvi (4-4) 

satisfies the KdV equation apart from a term of order rf. Thus it satisfies the whole 
hierarchy (still apart from terms of order rj 2 ). By linearization of the n th equation of the 
hierarchy (|3.1U|) . we find that 

d Tn vi = dxVnVx, (4.5) 

with 

V nVl = (d^ 71 - 1 + Ld x L n - 2 + ■■■ + v + C n vi (4.6) 

and 

dC 1 f x 

dl= d^ iVl) = ~ Vl + 2j dX ^ dxV ^- ( 47 ) 

This result is in particular valid for v \ = tp^ , but also for the following orders v\ = ip^ , 
with analogous notations, for any even p ^ 4. While the higher order time evolution 
of the part ip^ of the term of order e 4 coming from the right hand member H4 ((^2) 
of the linearized KdV equation is defined by the main order itself, the evolution of the 
part (p^ coming from the initial data is determined in an analogous way as the main 
term: <p^ satisfies the linearized KdV hierarchy ()4.5|) . 

Formal identities are related to this feature. They are found by the following reasoning. 
We assume as previously that v\ satisfies the homogeneous KdV equation. We have seen 
that v\ must satisfy equations (|4.5|) . The existence of v\ as a function of all the variables 
X, Ti,T2, . . . is not in doubt, because it is ensured by the existence of the Hirota r-func- 
tion [2j. Thus we have 



d Tn d Tj vi = d Tj d Tn vi 



(4.8) 
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for any integers n and j. We denote by {dT n Dj) the T n -partial derivative of T>j which de- 
pends on T n through v, and by dT n T>j the operator that applies successively T>j and dr n - 
With analogous notation, dxT^j = (dx"Dj) + T>jdx- With the use of equation (|4.5|) equa- 
tion (|4.8|) becomes 

((d Tn V 3 ) + VjOxVn) Vl = {(d Tj Vn) + V n d x V 3 ) Vl . (4.9) 

As in the case of the KdV hierarchy itself (nonlinearized), equation Q4.9JI is valid for any 
function v\. Thus the following identity holds formally: 

(dT n Vj) + Vjd x V n = {d Tj Vn) + V n d x V 3 . (4.10) 

This identity can also be written as 

[d Tj - d x Vj , d Tn - d x v n ] = o, (4.ii) 

where [M, N] = MN — NM denotes the commutator of the operators M and N, and 
the Schwartz conditions [d^ , dr n ] = are assumed. The commutator in equation (|4.11|) 
is easily computed and vanishes due to the Schwartz condition and to the formal iden- 
tity iniinj). 

Using the scaling definition for vq and T± , T2, ■ ■ ■ , we can write down the equations of 
the linearized KdV hierarchy for the functions ipp as 

+ (-4)™ 7 „9 € P n ^ 1 ) = 0, (4.12) 

with T> n defined by H4.6I) - ()4.7|) . In the case that n = 2 and p = 4 this equation has the 
explicit form 



5^72 
3 7 
5 /? 2 72 
6 7 2 



$dt<pP+2<p 2 (dt<p 2 )<pP]. (4.13) 



4.2 All secular terms vanish 

e fu 

if ip^ appears in H 6 (<p^,ip 2 ). Equation (|4.13|) yields 



The function cp^ is secularity producing [15 . The following term tp$ would thus be secular 



d T2 ^ ] = -729|rf ) +0 2 . (4.14) 

Together with equation Q4.2[) , this allows one to compute the linear part of the right hand 
member to obtain 

H 6 fa, P2) = -d^f ~ 72^4 2) + °2- (4.15) 

We admit that the only secularity producing terms due to (p^ are the linear ones and 
(2) 

recall that <p\ is not secularity producing. Thus, due to the equation of the linearized 



340 



H Leblond 



KdV hierarchy (|4.13jl . all secularity producing terms vanish from the right hand member 
Sg (ty?4, ^2) and ifQ is bounded. 

In order to justify the same property at any order, we have to compute the linear terms 
in ifp that appear in the expression of 3 n (</?2, • • • , The computation of all linear 

terms in these right hand members is detailed in the Appendix. The general expression 
of the linear part of 5 n is 

3 n fa, ... , (pn-2) 

£ 3 (K-)i>i s k) n / d Tj d\y v + TZ. (4.16) 

2jmj+fc+p=n+3 

The remainder 7£ contains nonlinear terms in </? p , as well as terms in (fj, j ^ p. From 
equation (|4.12|) of the linearized KdV hierarchy we see that 

d T] ^ = - lj df +l ^ + 2 . (4.17) 

Thus 

E n fa, ... , Vn-2) 

Sj-^i 2jmj+fc+p=n+3 

The remainder TZ contains nonlinear terms, terms in ipj, j 7^ p, and in the part (f p 
of (p p that comes from the right hand member 3 p (ip%, • • • > <Pp-2)- These latter terms 
are not secularity producing. Using the fact that, in the sum in (|4.18|) . m,j = for 
j > s = (n + 3 — p — 1) /2, and expression (|3.19j) of 7 j, we find that 3 n (992, • • • , { Pn-2) 
reduces to TZ' . This result is valid for each p, thus 3 n (</?2, • • • , { Pn-2) does not contain any 
secularity producing term and ip n is not secular. 



4.3 The higher order time derivatives 

On the other hand the expression of the solution 994 of the linearized KdV equation 
involves its spectral transform and an explicit computation of its T2-derivative seems quite 
impossible. However, it is needed for the explicit computation of the right hand member 
36 (^2, Wi)- d T2 ip^ is directly given by equation ()4.13|) . We now seek an analogous 

(2) 

expression of d T fa>. With the use of the definitions of this section, the linearized KdV 
equation (|2.42|) can be written as 

(8 T1 + 4 7 «%Di) tpf = 3 4 fa) ■ (4.19) 

We apply the operator (d T2 — 16 7 2<9gD 2 ) to both sides of this equation and make use of 
identity ()4.11j) . to obtain 

(8 T1 + 4 7 %£>i) (8 T2 - 16 72 a e P 2 ) ipf = (d T2 - 16 72 %X> 2 ) 3 4 fa) ■ (4.20) 
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(<9 T1 + Ajd^T>i) is the operator on the left hand side of the linearized KdV equation that 

(2) 

admits a unique solution for a given initial data. tp\ is defined in such a way that it 
vanishes at t\ = for any £ and t<i- Thus its £ and T2 derivatives also vanish and 



(d T2 - 167 2 0fE> 2 ) <p K i 



(2) 



71=0 



0. 



We use the following notation to denote this solution: 

(d T2 - 16 72 %2? 2 ) <pf = (d T1 + 4 7 c Pi)- 1 p T2 - 16 72 ^r> 2 ) ~ 4 (^ 2 )] . 



(4.21) 



(4.22) 



d T2 E,4 (yj 2 ) is easily computed while the expression of H 4 is explicitly known through the 
perturbative expansion and d T2 (p2 is determined by the second equation of the hierar- 
chy ipTTBJ) . It is 



a - i \ ic ^2*4 (^2) fl ~2 
(9 T2 = 4 (V?2j = -I672 j <%£ <^2, 

a<^ 2 



(4.23) 



where dH 4 ((p 2 )/dif2 is the differential operator obtained by linearization of H 4 (<p 2 ) about 
(/?2 • 9 T2 (f4 is thus given by the formula 



d T2 v4 2) = 1672 ( d(D 2V >f - (d T1 + ^d^y 1 



dZ 4 (ip 2 ) a r2 



dif2 



(4.24) 



(2) - 

More generally tp p is the solution of the linearized KdV equation with the right hand 



member E p (<^ 2 , • • • , fp-2) and zero initial data and 

(d T1 + Ar,d(D x ) (d Tn - {-A) n ln d(D n ) <pV> 

= (9 Tn - (-4) n 7 n a c P n ) E p (ip 2 , (pp- 2 ) ■ 



(4.25) 



Because, for n 7^ 1, (d Tn — {—4) n ^ n d(V n ) does not contain the partial derivative d. 
{d Tn -{-A) n ln d t V n )4 



71 ! 



(2) 



n=0 



0. Using the above notations we obtain the following 



expression of the higher order time derivative: 

,<Pp-2) 



"p-2 

£ 

.21=2 



d<p. 



21 



d Tn V2i - (-4) n 7 n ^X» n H p (<p 2 , 



(4.26) 



For p = 4 the time derivatives in the right hand member of equation (|4.2fiJI reduce to 
Qr n f2 = —(—^) n 7nd^C n ip 2 , according to equation (|3.18l) . For larger values of p, it involves 
d Tn (f2i with 2/^4 that divides into d Tn p 2 i = 0^^' + d Tn ip^ , where c^^j * s gi ven by 
the linearized KdV hierarchy (|4.12|) . and d Tn ip 2l by the same equation (|4.26j) . in a recurrent 
way. 

We have proved that the introduction of the KdV hierarchy removes all unbounded or 
secular solutions from the perturbative expansion and given the expression of the higher 
order time derivatives of all terms in the perturbative expansion. 
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5 Conclusion 

We have studied the higher order terms in the perturbative expansion that describes 
KdV solitons in a rather general frame, including electromagnetic wave propagation in 
ferromagnetic media. Using various mathematical techniques, we have been able to write 
down the equations satisfied by the quantities of any order in this expansion. In every case 
it is the linearized KdV equation, with some right hand member. We have summarized 
the known results about this equation. Unbounded or secular solutions can be removed 
by suppressing linear terms in the right hand member of the equations. This is done by 
imposing that the main term satisfy all equations of the KdV hierarchy. There exist scaling 
coefficients for the higher order time variables in the hierarchy, which are determined by 
the requirement that the linear terms in the right hand member of the linearized KdV 
equation vanish. They are computed using explicit recurrence formulas. Furthermore the 
behaviour of the higher order terms relative to the higher order time variables was not 
known. These terms are divided into two parts: one part comes from the initial data, and 
its higher order time evolution is described by the linearized KdV hierarchy. The other 
part comes from the right hand member of the linearized KdV equation that describes its 
evolution relative to the main time scale. It is completely determined, including its higher 
order time evolution, by this latter equation. 

A quasi-explicit formula for the higher order time derivatives of this second part of the 
corrective terms has also been given. It involves also the expression of the linearized KdV 
hierarchy. The main result is that none of the higher order terms produces any secularity. 
Thus the existence of the expansion is established up to any order. 
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A Appendix. Computation of the linear terms 
in the right hand member of the equations 
of the perturbative expansion 

The right hand member 3/ (tp 2 , <f3, ■ ■ ■ , fi-x) of the linearized KdV equation (|2,42|) is 
defined by 

a Ei ((p 2 , ip 3 , . . . , ^_i) = ni<Sj (i+3) + a (d T1 ipi + j3d^{ip 2 ipi) + 7<9f ¥>0 (A.l) 
for each I ^ 3, with 

n-2 

£q^ n) = -Eu n - {A- V)d^u n ^ 1 - )d T .u n - 2 j-i + B(u p , u n _ p ), (A.2) 

3>1 P=2 
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so that the n th equation (|2.34j) - (|2.3fij) of the perturbative scheme is <S?( n ) = (&/(«) is 
the difference of the right hand member minus the left hand member). We write O2 for 
any polynomial expression in the ipj and their derivatives and primitives that does not 
contain any linear term. This is analogous to the usual Landau notation 0(</?|), except 
that we do not assume that the ipj are small in any way. Note that the use of the KdV 
equation satisfied by ip2, and so on could change the linear terms in H; (cp2, ^3, • • • , Pi-i)- 
Thus in order to define the linear part of 3; (ip2, </?3, • • • , Vz-i) m a unique way, we assume 
that no use of these properties has been done. We made the standard hypothesis that, 
for zero initial conditions at these orders, all ifj with an odd value of j are zero, but 
it is more convenient for the present formal computation to forget this feature. Formal 
first order terms, vector u\ and function ipi, satisfying equations ()2.34[) and (|2.35|) are 
introduced, although it can be shown that they are necessary zero. This is consistent with 
the requirement that no use of the equations has been made: we compute formally the 
right hand member without solving the equations in any way. 
We write an a priori formula for the linear part of u n : 

u n = u(m)HU c) T .J dlw+Oi, (A.3) 

m=((mj)j^i ,k,l) 3*^1 
m,j,k^zO,l^l, d(m)=n 

where the components of the vector coefficients u(rh) have to be determined. The homo- 
geneity properties of u n are described by the "degree" d(fh), defined by 

d {{mj) m , k, I) = 2 i m 3 j +k + l (A.4) 

No explicit dependency on n has to be written down because n = d{rh). The substitution of 
formula ()A.3|) into the recurrence formulas ()2.35j) and l)2.36j) for u n yields a new recurrence 
formula that allows one to compute these coefficients. This shows by induction that 
formula (|A.3|) is valid for any value of the integer n and that the values found for the 
coefficients are valid. These recurrence formulas are 

For all 1^1: u ((0), 0, 1) = a (A.5) 

For all k and I ^ 1: u ((0), k, I) = F^u^O), k - 1, 1)). (A.6) 

For all (rrij)j^i / (0) and / ^ 1, 

u ({ mj ) m ,0,l) = ^F 3 (n ((m,- -5ij)j>i,0,0). (A.7) 

For all (mj)^i ^ (0), k, I ^ 0, 

u ((mj)j^,k,l) = F^u^rrij)^!^- 1,1)) 

+ ^2 [F 2 {u ((rrij - 5 itj ) m , k - 1, 1)) + F 3 (u(( mj - 5 itj ) m , k, I))] . (A.8) 

i>l 
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8ij is the Kronecker symbol and the operators Fj are defined by 

F x =F 2 (A-V), (A.9) 
F 2 = — [l — (A — V^QM - V)} E- l Q , (A.10) 
F 3 = -(A-F)- 1 Q 1 . (A.ll) 

Formulas (|A,5|) to (|A,8|) define u ((nij)j^i,k,l)). We see that this quantity does not 
depend upon I. We write more simply u k)) in the following. For some particular 

terms explicit expressions can be given. It is straightforwardly seen from equations ()A.6|) 
and (|A.5|) that, for all k, 

u((0),k) = Fi fc a . (A.12) 
In the particular case that k = 1 this is 

u ((0),1)) = ai +a[. (A.13) 
Equations (|A.7|) and (|A.5|) yield in a similar way 

u(K0 3 >i,0) = F 3 ^i m ^o. (A.14) 

The definition of ao shows that Qiao = 0. Thus F^ao is zero and u((mj)j^i, 0) = for all 
( m i)j>l- I n the same way, due to Qoao = 0, i^ao = 0. 

The formulas obtained are substituted into the definition (|A.2|) of to give 



d(fh)=n— 1 



£ ni«(m)n(/ M d k +\ l + <D 2 . (A.15) 



(mj)j>i,fc^0,i^l ' 



d(m)=n— 2i— 1 

Thus the right hand member 3 n is 

(^2, f3, • • • ¥>n-l) 



2 Sam^.^n^^.J %>i + 02, (A.16) 



l<j <n-l 
Sj^i 2j"rry+fc+i=n+3 



with 



a S ((mjj^i, fc) = -IIi (A - V)Q Q u ((mj)j^i, k - 1) 

- n x « ((m 3 - - Sij) j>lt k-1), (A.17) 



■t>i 
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where u((mj)j^i,k) is defined by the recurrence formulas (|A.5|) to (|A.8|) . Note the re- 
markable feature that the coefficient E (('mj)j^i, k) does not depend on n. 

The term in which d Tj ipi appears with the largest value of the index j, in a given E n , 
can be computed explicitly. Because n is necessary even, we write n = 2p. Due to the 
homogeneity properties of E n (ipi, . . .) we see that the term sought is proportional to d Tp ip2- 
The corresponding coefficient is E ((Sj tP ), 1). By direct application of the previous formulas, 
we find that 



S((«5 i>p ),l) = -1. (A.18) 

The coefficient 72 that defines the second order time scale ti is computed using the ex- 
pansion (|3.8|) and the above recurrence formulas. It is 

72 a = Hi [-7F1 + (A - V)Q (Fi 3 - 7 [F X F 3 + F s Fi + F 2 })} F iaQ . (A.19) 
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